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Abstract 

Functional Renormalization Group Equations constitute a powerful tool to encode the perturba- 
tive and non-perturbative properties of a physical system. We present an algorithm to system- 
atically compute the expansion of such flow equations in a given background quantity specified 
by the approximation scheme. The method is based on off-diagonal heat-kernel techniques and 
can be implemented on a computer algebra system, opening access to complex computations in, 
e.g., Gravity or Yang-Mills theory. In a first illustrative example, we re-derive the gravitational 
/3-functions of the Einstein-Hilbert truncation, demonstrating their background- independence. 
As an additional result, the heat- kernel coefficients for transverse vectors and transverse-traceless 
symmetric matrices are computed to second order in the curvature. 
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1 Introduction 



1.1 The Functional Renormalization Group 

Functional Renormalization Group Equations (FRGEs) constitute a universal and highly flexible 
tool for unlocking non-perturbative information in a plethora of settings ranging from condensed 
matter physics, to statistical and elementary particle physics, to quantum gravity jll EJ [3l HJ 
[5]. They provide an explicit realization of Wilson's ideas of renormalization [6j, allowing us to 
integrate out quantum fluctuations in a system "shell-by-shell" in a much more manageable way 
than, e.g., directly manipulating the path-integral. 

One of the main practical virtues of the FRGEs is that general non-perturbative properties 
of a theory, like phase transitions or fixed points of the renormalization group (RG) flow, can 
be found via relatively simple approximations. This feature makes FRGEs a primary tool in 
the search for non-Gaussian Fixed Points (NGFPs) of the RG flow in Quantum Field Theories, 
which provide a non-trivial and predictive generalization of the notion of perturbative renormal- 
izability within the framework of Asymptotic Safety [71 18]. The existence of such NGFPs has 
been well-established in the context of, e.g., 0(A^)-sigma-models [9] and has recently also been 
discussed in the context of standard model physics related to Yukawa-systems [ini HI [IS] and 
the Higgs sector [T3], where it has lead to new physics scenarios beyond the ones allowed within 
perturbative renormalization. A NGFP has also been conjectured to control the UV behavior of 
four-dimensional gravity [T6| (see also |17l 118]). rendering the theory non-perturbatively renor- 
malizable or, equivalently, asymptotically safe. This conjecture is motivated by the fixed-point 
structure of gravity in d = 2 + e dimensions [141 US] 3.nd evidence in its favor has been established 
in a series of works (see [191 120] |2T1 122] for reviews and detailed references). 

While simple approximations of the exact RG flow, usually consisting of truncating the flow 
to a handful of running couplings, can already yield qualitative information on the fixed-point 
structure of a given theory, establishing the consistency of the approximation and answering ques- 
tions concerning, e.g., critical exponents, universality classes and the number of relevant couplings 
appearing in the theory requires more sophistication. For the case of gravity, in particular, the 
best truncations to date [231 ISl (SHI [261 1211 [28] have mostly been restricted to fluctuations around 
a spherical background, and are insufficient to properly answer these questions. In this light, it is 
highly desirable to have a computational tool which can track the RG flow on a sufficiently large 
subspace of all running coupling constants. 

The main purpose of this paper is to present such an algorithm, which we call the "universal 
RG machine", that allows us to compute approximations of the exact FRGE in a flexible and 
systematic way. Its key element is the algorithmic reduction of the (usually complicated) opera- 
tor traces appearing in the flow equation to the off-diagonal heat-kernel expansion of a standard 
Laplace-operator on a curved manifold [23 [30] (see also [311 [32] for current results and further 
references). This procedure is general and does not rely on the choice of a particular background, 
which is usually otherwise evoked to simplify computations. Furthermore, the bookkeeping re- 
quired in a practical computation can be easily handled by a computer algebra program, opening 
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the door to computations which until now have been out of reach due to technical limitations. 
The algorithm will be presented here having gravity as the underlying theory, but we stress that 
the construction is readily adapted to other gauge-theories. 

The rest of the paper is organized as follows. In the next subsection, we give a short introduc- 
tion to the Wetterich equation [33], which provides the foundation of our construction. Sect. [2] 
contains the blueprint of the universal RG machine and constitutes the main result of our work. 
As an illustration, we apply our construction to the Einstein-Hilbert truncation [5 ^ [55 1 [55 1 [57 1 [55] 
in Sect. [3l before closing with an outlook on future research perspectives in Sect. HI Proofs of the 
more technical details entering the construction have been relegated to five appendices, which 
cover the commutator relations, transverse-traceless decomposition of the fluctuation fields, the 
derivation of the heat-kernel coefficients for transverse vectors and transverse-traceless matri- 
ces, and general considerations on the structure of the flow equation regarding gauge- flxing and 
truncation independent contributions. 



1.2 The Wetterich Equation 

The starting point of the universal RG machine is the FRGE for the effective average action 
[33], which schematically takes the form 



dtrk[<PA] = iSTr 



+ TZk) ' dtTZk 



5, = . A. (1.1) 



Its construction for gauge-theories [59j and gravity [51] employs the background field formalism, 
which ensures the background covariance of the flow equation and implies that Tk[(t)-, depends on 
two arguments, the background fields (j) and the fluctuations around this background, collectively 

(2) 

denoted by 4>. The Hessian T^ is the second variation of with respect to the fluctuation fields 
<pi at fixed background (j) and is, in general, matrix-valued in field space, 

1 1 5^T 



r(2) 
^ k 



\/5 Vg HiHj ' 



where g is the determinant of the background space-time metric (which equals unity for flat 
space) and [j] takes the values 0,1 for bosonic and fermionic flelds respectively. In the FRGE, 
STr is a generalized functional trace which includes a minus sign for ghosts and fermions, and 
a factor of two for complex flelds, while TZkip^) is a matrix valued IR cutoff which provides a 
fc-dependent mass-term for fluctuations with momenta < k'^. The cutoff satisfles TZ^ oc k'^ 
for p2 <^ k'^ and vanishes for high-momentum modes, TZk as ^ k^. The factor dtTZk 
ensures that the r.h.s. of the flow equation is flnite and peaked at momenta ^ k'^, so that an 
additional UV regulator becomes dispensable. Eq. (jl.ip describes the dependence ofTk[4>,4'] on 
a renormalization group scale k, realizing the Wilsonian picture of renormalization by integrating 
out quantum fluctuations with momenta k'^ < < A. Its solutions interpolate between the 
ordinary effective action T[i^\ = Tk=Q[4' = 0, i;^] and an initial action Fa at the UV scale, which in 
the limit A — )• oo essentially reduces to the bare action (see ^40j for more details). 
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The main shortcoming of the FRGE is that it cannot be solved exactly. In order to extract 
physics from it, one thus has to resort to approximations. In the one-loop approximation, where 
Ffc on the r.h.s. of (jl.ip is replaced by the /c-independent bare action, perturbative results are 
reproduced. Going beyond perturbation theory, a common strategy relies on expanding rfc[(/), 
on a basis of local invariants and truncating the series, thereby limiting the form of the effective 
average action to a finite number of invariants In[4', <P], 

N 

rk[^J]^^Un{k)In[^J]. (1.3) 
n=l 

Substituting such ansatz into and projecting the exact RG flow onto the subspace spanned 
by the In[4>,4'] then provides the /3-functions for the fe-dependent couplings Un{k). As this ap- 
proximation scheme does not rely on an expansion in a small parameter, the resulting /3-functions 
typically capture contributions from all loop-orders and constitute non-perturbative results. Pos- 
sible ordering principles for organizing the expansion include, e.g., the derivative expansion, 
whereby interactions with increasing number of derivatives are systematically included in the 
ansatz, or a vertex expansion, in which (|1.3p is truncated at a certain order of the fluctuation 
field (/.. 

Depending on the complexity of the ansatz p.3p . projecting the right-hand side of the flow 
equation onto the In[4'^ 'P] can become quite involved, since the matrix- valued operator (F^ +Tlk) 
is usually of non-minimal form, which poses a severe complication for the evaluation of the 
functional trace in (jl.ip . This obstacle is frequently circumvented by working with a particular 
background configuration, with the drawback that some In[4>,4>] become indistinguishable, thus 
potentially hiding important features of the theory under investigation. In the gravitational 
context, for example, it is common to restrict the background metric to be maximally symmetric, 
i.e., a sphere. In such a case the only terms that can be recognized are integral powers of the 
Ricci scalar, and it is not possible to distinguish other important invariants, such as the terms 
/ = / d'^x^R^"'Rf,^ or / = / d'^x^RD'^R, which can affect the counting of the degrees of 
freedom of the theory. In order to distinguish such invariants, we need a generic background, but 
this leads to serious technical complications due to the operator structure appearing in the second 
variation of the action. It is thus clear that a new computational scheme, enabling us to deal 
with the general operator structures appearing in the FRGE without resorting to a simplifying 
choice of background, is highly desirable. In the next section we will propose such an algorithm 
in form of the universal RG machine. 

2 The universal renormalization group machine 

We will start by discussing the universal RG machine from a general perspective, before giving 
an explicit example of its implementation in the case of gravity in Sect. [31 The algorithm that 
we propose consists of the following steps: 
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1. expanding to quadratic order in the fluctuation fields (pi; 

(2) 

2. simplifying the operator structure of ' by performing the transversal 



decompositions (j2.4p and (|2.6p of the vector and tensor fluctuations, 
respectively; 

(2) 

3. decomposing into its kinetic part IfC and interactions M according 
to dsn; 

5. inverting [F^ +Tlk] via the perturbative expansion in the M's, eq. (|2.12p : 

6. inserting projection operators in the traces; 

7. bringing the trace arguments into the canonical form (j2.15p by applying 
commutator relations for covariant derivatives; 

8. evaluating traces without interactions via standard heat-kernel expan- 
sions for transverse fields; 

9. replacing derivatives appearing in M by the ff-tensors (|2.23p constructed 
from the off-diagonal heat-kernel; 

10. reading off /3-functions by matching background structures on both sides 
of (ini) ; 

Notably, each of these steps can be handled by computer algebra software. This is of crucial 
advantage when dealing with more complex expansions (truncations), since these require an 
increasing amount of bookkeeping when collecting the various contributions to the RG flow. We 
now proceed by discussing the steps indicated above in more detail. 

(2) 

2.1 Simplifying the operator structure of 



k 



Schematically, f!^^ can be written as 



r(2) 
k 



i{A)6'^h + n'^iDf,) + M'^iR,D^) , (2.1) 



kinetic terms uncontracted derivatives background curvature 

where i,j denote the indices in field space, Ij denotes the unit of the corresponding field space, 
and A = —g^''^D^Djj is the Laplace-operator. The first piece, ICj(A), is the diagonal part of 
the inverse propagator of the i-th field, containing only Laplacians and eventually cosmological 
and mass terms. The last two pieces group the off-diagonal operators, with M*-^ including all 
interactions that involve at least one power of the object controlling the expansion of the flow 
equation (11. Sp . and with B*-^ containing all the remaining terms with derivatives but no such 
expansion object. For a single metric truncation of the flow equation, the expansion field is 
naturally given by the curvature tensor of the background metric. The formalism, however, 
allows for an expansion in any background structure like the background ghost fields or, in 
the case of a vertex expansion, the fluctuation fields </>. Explicit examples of typical off-diagonal 
matrix entries are 

D = (1 - 0)0^,0" , B = D'^D'^DaDp , (2.2) 

and 

M = R'^'D^D^ , M = Df'cu . (2.3) 
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Our goal is to evaluate the functional trace on the r.h.s. of (jl.ip by means of off-diagonal heat- 
kernel techniques. The first step to this aim is to simplify the operator structure of (j2.ip so 
that terms of the D-type do not occur. The reason for this is that operators of this kind would 
contribute with arbitrary powers, spoiling the perturbative inversion of the modified propagator 

(2) 1 

-|- TZk which we will subsequently employ. 
One way to remove these operators, often used in connection to ordinary perturbation theory, 
is to choose a convenient gauge- fixing condition. Considering the first example in (|2.2|) . for 
instance, the corresponding B-term could be removed by choosing the value a = 1 for the gauge 
parameter. Following this route is problematic for two reasons, however. Firstly, it restricts the 
gauge-sector to a very particular choice, which has to be carefully adapted to the truncation 
to be studied. Secondly, in most cases there will not be enough freedom in the gauge-fixing 
function to remove all B-type terms. The latter is the case for, e.g., higher-derivative gravity, 
where conveniently-chosen gauges can remove all four-derivative B-terms, but lower-derivative 
terms still remain [41j. 

A systematic way to remove these terms is the transverse-traceless- (TT) (or York-) decom- 
position of the fluctuation flelds on a generic gravitational background. For vector fluctuations 
this decomposition is given by 

ct>^, = cl)l + D^iq , (2.4) 

subject to the differential constraint 

D>^<t>l = ^. (2.5) 
Analogously, metric fluctuations are split according to the "minimal" TT-decompositior0 

V = + D^iu + D^^f, - Ig^.uD'^ia + \g^,uh , (2.6) 

with the component fields satisfying 

D%u = ^. 9'"'hl, = 0, g^''h^u = h. (2.7) 

At the price of introducing additional fields, the transverse constraints on the fiuctuation fields 
guarantee that any B-type term can be converted into M-contributions. This follows from the 
fact that at least one covariant derivative contained in D has to be contracted with a transverse 
fluctuation fleld. Using standard commutation relations, this derivative can then be moved to the 
very left or very right of D so that it acts on the fiuctuation in such a way that the resulting term 
vanishes. In the process, one picks up additional terms which contain background curvatures and 
are therefore of M-type by definition. As a result, after performing the TT-decomposition of the 

(2) 

metric fiuctuations and rearranging the resulting operators, we end up with no B-terms in T^ . 



^The "standard" TT-decomposition further decomposes the vector according to (|2.4|) . as is explicitly done 
in App. [E] 
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2.2 Constructing the inverse of 



(2) 

Having obtained the inverse propagator T) in simplified form, the universal RG machine proceeds 



by expanding the latter in 



with the construction of the regulated propagator 
the background structure controlling the truncation, as, e.g. the background curvature. The 
central idea here is to treat the off-diagonal, M-type operators appearing in the inverse regulated 
propagator as perturbations on the diagonal, Laplace-type operators in K. The presence of the 
expansion field in the M-terms allows us to perturbatively invert + TZ^ in M along the same 
expansion scheme as that of the effective average action, thereby accounting for all contributions 
to any given truncation. 



2.2.1 Adapting the infrared cutoff 



(2) 

Following standard FRGE practice, we first modify the inverse propagator F^ by a cutoff TZk- 
With the structure of the inverse propagator split into K and M, it is natural to tailor the cutoff 
to the K-term alone. Following this scheme, the Laplacians appearing in K are dressed up by the 
momentum-dependent IR regulator, according to the rule 



A^Pk = A + Rk{A) 



(2i 



where -Rfc(A) denotes the basic (scalar) cutoff function, which interpolates monotonously between 
Rk{z) oc k'^^z <^ k'^ and Rk{z) =0,z ^ k'^ . As a consequence, the matrix-valued kinetic terms 
in the modified inverse propagator take the form 



]fC(A) ^ P(A) = ]IC(A) + Mfe(A) = K{Pk) 



(2.9) 



where P is obtained from K by applying the replacement rule (j2.8p . This construction guarantees 
that the IR-regulator IZ^ is diagonal in field space. Note that K, and therefore also P and M,fc, will 
in general contain A;-dependent coupling constants such as wave-function renormalization factors. 

It is useful to compare our construction with cutoff scheme classification of [38]. Since ()2.8p 
regulates the Laplacian only, this scheme falls into the class of Type I cutoffs. However, since 
only Laplace-operators appearing in IK are regulated (and not those that could occur in M) this 
is not the standard Type I cutoff, but a variation of it which we will call "Type Ic" scheme in 
the following. We note that the Type I and Ic cutoffs agree for truncations containing only two 
space-time derivatives, while they will generically differ for higher-derivative truncations. 



2.2.2 Inverting the regulated inverse propagator 

Having completed the cutoff implementation, the inverse regulated propagator becomes a block- 
matrix in field space 



F(2) + Uk 



P^(A) 1a + Ma 



Pb(A) 1b + Mb 



(2.10) 
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While we restricted ourselves to the case of two fields here, as this is what is typically encountered 
in practical computations, the generalization to the case of more than two interacting fields is 
straightforward. 

In order to obtain the modified propagator from (j2.10p . we first formally apply the exact 
inversion formula for block-matrices, 



" A 


B ' 


-1 


C 


D 





{A-BD-^C) 
-D-^C [A - BD-^C) 



-A-^B {D 
(D-CA 



CA-^B) 



(2.11) 



Building on the split between the kinetic terms P and the interaction terms M, each block-matrix 
entry above can be constructed perturbatively in mJ§ This computation may be simplified further 
by observing that the cutoff operator TZi~ and thus dtTZk are diagonal in field space. Thus, it 
suffices to consider the diagonal elements in (|2.1ip . Substituting for (j2.10p . the perturbative 
inversions in M for the diagonal elements read 

-1 



r(2) + TZk 



.(2) 



AA 
-1 

BB 



-Ma — + —Ma — Ma — + —'. 

Fa Fa ^ Pa Pa Pa Pa Pa 



1 



it ^ + 0(m3) 



1 



1 



■Ub 



1 



+ 



1 



1 1 

P^^P^ 



1 



1 



+ 0{I 



(2.12) 



where the Mj are understood as suitably projected onto the subspace spanned by the fluctuation 
fields. As M contains the background quantity organizing the expansion scheme of the flow equa- 
tion, we are ensured that only a finite number of terms appearing in these expansions contribute 
to a given truncation. 



2.3 Evaluating the traces via off-diagonal heat-kernel expansion 



Substituting the perturbative expansion of 



into the r.h.s. of the flow equation 



results in a series of individual operator traces on a space of fields satisfying the differential 
constraints imposed by the TT-decomposition. We will now discuss the evaluation of these 
traces, distinguishing them into two different classes according to whether or not they contain 
background vertices M. 



2.3.1 Traces without non-minimal operator insertions 

The traces without background vertices typically assume the simple form 

(2.13) 



^We stress that this perturbative inversion should not be confused with standard perturbation theory, where the 
power series is controlled by a small coupling constant. Here, the perturbative series is controlled by a background 
quantity, the expansion field, which organizes the truncation of the effective action. Consequently, this approach 
still allows us to work "non-perturbatively" in the coupling constants. 



^kin 



TV, 
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U 


i 


o 
z 


i i 


r)rp 

z i 


Cl 


1 
i 


/I 


1 n 

iU 


Q 
O 





C2 


1 
6 


2 
3 


5 
3 


1 

4 


5 
6 


C3 


1 

72 


1 

18 


5 

36 


24 ^ 2)^E 


137 1 N 
216 2)^E 


C4 


1 
180 


1 

45 


1 

18 


1 2 
40 


17 2N 
108 


C5 


1 
180 


11 
180 


4 
9 


15 ^ 2ye 


5 _|_ Af 
18 2yE 



Table 1: Coefficients appearing in the heat-kernel expansion (j2.14p for d = 4. The coefficients 
a for the IT and 2T-fields are novel and computed in Appendix O Here, the terms including 
the Euler characteristic xe are due to zero modes of the decompositions, and is the sum of 
conformal and non-conformal Killing-vectors of the background manifold. 



where s indicates the "spin" of the fluctuation field and s = 0, 1,2,1T,2T denotes the trace 
with respect to a scalar, vector, symmetric tensor, or a vector and symmetric tensor satisfying 
transverse-traceless conditions, respectively. As this class of traces contains only Laplace opera- 
tors in their arguments, they can be evaluated via standard heat-kernel methods by expressing 
the operator occurrence as a Laplace-transform and then evaluating a standard heat trace of the 
form 

Tr, [e-^^] = j^^^ J d'^x^ [ci + C2R + OiR" + c^R^.R^'' + c^R^.afsR''"'"'] ■ (2.14) 

For fields without differential constraints, these coefficients are well-known and can be found 
in, e.g., \i2\ I43j. For fields satisfying differential constraints, IT and 2T coefficients are known 
for the special classes of maximally symmetric backgrounds [36j or Lichnerowicz-Laplacians on 
Einstein-spaces [2 7) only. In Appendix [Ul we have computed these heat kernel coefficients for the 
case of a general, purely gravitational, background. The result of this computation is summarized 
in Table [T] for d = A. 

2.3.2 Traces containing non-minimal operator insertions 

The second kind of contributions to the flow equation arises from traces that contain background 
vertices M with derivative operators of non-Laplacian form. In order to evaluate these, we first 
employ commutation relations for these covariant derivatives to arrange the trace arguments into 
"standard form" , so that all Laplacians are collected in a single function and the product of the 
background vertices forms a single operator insertion O. The trace can then be written as 

S^^'^^^ = Tt,[W{^)0] . (2.15) 

Denoting the Laplace transform of H^(x) by W{s), this trace can be re-expressed as 

l-OO 

Tii[W{/\)0]= dsW{s){xi\(i-'^0\xi) , (2.16) 
^0 
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where \xi) form a basis of eigenfunctions of the Laplace-operator in position space, carrying 
appropriate Lorentz indices. Inserting a complete set of states, ()2.16p can be evaluated using 
off-diagonal heat-kernel methods |32] . Explicitly, 

{x,\Oe-'^\xi)={xi\0\x[){x',\e-'^\xi) = j d^x^tn [OH{s,x,x')]^^^, , (2.17) 
where H{x,x';s) is the heat-kernel at non-coincident points 

/ °° 

H{x, x'; s) := (x'|e-^^|x) = (Atts)-'^/^ e""^ ^ s"A2„(x, x') . (2.18) 

n=0 

Here, a{x, x') denotes half the squared geodesic distance between x and x' and A2n{x., x'), n G N, 
are the off-diagonal heat-kernel coefficients. The latter are purely geometrical two-point objects, 
formally independent of the space-time dimension d, H and explicit values for them can be found 
in, e.g., |44] and references therein. 

In practice, the explicit evaluation of the trace (j2.17p requires that the covariant deriva- 
tives contained in O and acting on H{x,x';s) be totally symmetrized in their indices. This can 
be achieved by successively writing non-symmetric combinations of the derivatives as a sum of 
symmetric and antisymmetric pieces and expressing the latter in terms of curvatures using the 
commutation rule ()A.ip . The coincidence limit of the tensors arising after acting with the sym- 
metrized derivatives on H{x,x';s) can be computed making use of the following properties of 
a{x,x') and A2nix,x'). Denoting the coincidence limit x — >• x' by an overline and using round 
brackets for symmetrization of indices with unit strength, (a/3) = ^(a/3 + /3a), etc, we have 



2 



f, = a, a{x,x') = a{x,x').f, = 0, a{x,x').^^ai...a„) = ^ , n > 2 . (2.19) 



The only non-trivial contribution from a{x,x') arises when exactly two covariant derivatives act 
symmetrically on it, 

(^{x,x').(^^^) = gf„y{x) . (2.20) 

Secondly, up to terms of 0{R'^), the contributions from derivatives acting on the off-diagonal 
heat-kernel coefficients A2nix,x') become 

1 ^ 1 



Ao{x,x') = l, Ao(x,x');^ = 0, Aoix,x'),^^^) = -Rf,^, A2ix,x') = -R. (2.21) 
Using these relations and defining 



Ha,...a,„ ■■= " " " ^a^^) H{x, x'; s) , (2.22) 

we arrive at the following expansions, valid up to terms of 0{R^) and given here for later reference, 

Har...a2n = (4^g^d/2 { ("^s)"" (^aiaa " " " 5a{2„_i)a2„ + " 1) perm.) [l + ifii?) 

+ g ( — 2s) *■ ^ (^9aia2 ■ ■ ■ 5'«2n-3"2n-2-^a2n-l«2n + ( 2"{n-l)l ~ Pemi.^ ^ 



(2.23) 



'For some early work on the spin-dependence of the An{x,x') see also [301 , 
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Explicitly, 

Hal3 



1 



(47rs)°'/2 
1 



(2.24) 



For a given truncation, the types of i/-tensors required can be determined by observing that each 
power of the curvature contained in M can be contracted with at most two covariant derivatives. 
To a given order n in the curvature, their total occurring number is therefore given by 2n for 
a scalar trace, whereas for traces over (unconstrained) vector fields and symmetric tensors, the 
number of required derivatives increases by 2 and 4 respectively, since they carry open indices. 
This number is, of course, modified if the background structure used in the expansion carries a 
different set of indices. 

Having the i?-tensors at our disposal, it is now easy to systematically evaluate the perturbed 
operator traces (j2.16p . Expanding 



(2.25) 



k=0 



with totally symmetric matrices M°i'"°2fe^ we have 



ds W{s)Tr, 



-sA 



•«2fc ]J 



k=0 



d'^x^ / dsW{s)tvi 



D. 



ai...a2k 



1 r°° 



.k=0 
I ■ 



(2.26) 



Here, Ii indicates the interaction monomials contained in the truncation ansatz, i.e., I2 
J d'^x^RJ-s = f d'^x^R^, etc. 

The final form of the result is most conveniently written in terms of the Mellin-transforms 



Qn[W] 



dss-'^Wis), 



which we can re-express in terms of the original function W as 



Qn[W] 



1 



Tin) 



dzz^-'W{z), Qq[W] = W{Q). 



(2.27) 



(2.28) 



In view of practical computations, we furthermore note that there is a degeneracy between Q- 
functionals if their argument contains powers of the Laplacian as a factor 



(2.29) 
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For a given truncation, the r.h.s. of the flow equation (jl.ip is then obtained by adding the 
relevant contributions of all the individual traces (j2.26p and (j2.13p . The desired /3-functions 
governing the scale-dependence of the coupling constants are then read off as the coefficients 
multiplying the interaction monomials spanning a given truncation. 

We close this section by highlighting the central virtues of the algorithm outlined above. 
Firstly, it is completely algebraic. No numerical integrations are needed in order to determine 
the /3-functions. Secondly, the shape of the cutoff function TZi^ is left arbritray. Thirdly, and most 
importantly, at no point in its construction does it make any reference to a particular choice of 
background. This is one of its main improvements compared to other computations. All previous 
explorations of the gravitational theory space, for instance, relied heavily on the fact that the 
background metric satisfies certain properties, as e.g., being the maximally symmetric metric on 
the d-sphere, or an Einstein metric. For the systematic construction of the derivative expansion 
in gravity theories, relaxing this technical limitation is of central importance. 

3 Example: The Einstein-Hilbert truncation 

We now illustrate the working of the universal RG machine within the simplest gravitational 
setting, the single-metric Einstein-Hilbert truncation in d = 4 [SH [351 ESI EZl [38]. In this case, 
the expansion scheme is controlled by the curvature tensor of the background metric and the 
truncation retains all terms up to linear order in R. In particular, this implies that terms which 
contain a covariant derivative of a curvature scalar (e.g. D^jRap^u, etc.) do not contribute 
to the truncation and can safely be neglected. Besides demonstrating that the new algorithm 
leads to reliable results, the main emphasis of this section is on demonstrating the background- 
independence of the resulting /3-functions0 Instead of working with a specific background metric, 
our construction will leave g^y unspecified and the only technical assumption is that 5^1^ is a 
(Euclidean) metric on a compact, closed, and complete manifold, which guarantees that the 
minimal TT-decomposition and heat-kernel expansion are well defined. 

Splitting the averaged metric according to g^j^ = g^u + /i^j^, the Einstein-Hilbert truncation 
corresponds to the following ansatz for the effective average action 



and the gauge fixing and ghost actions are defined in ()E.2p and ()E.13P respectively (see Appendix 
[E] for the derivation of these truncation- independent terms). For convenience, we combine the 
scale-dependent dimensionful Newton's constant Gk and cosmological constant A^ into 



Tk[K C, C- -g\ = rri5 + h] + S^'[h- g] + S^'^-'^h, C, C; g] + S' 
where the gravitational part r|^'^^[(jr] is 



laux 



(3.1) 




(3.2) 



Afe 



1 



(3.3) 



Mo = 



ui = - 



^See also [38| for a related discussion. 
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and introduce the dimensionless couplings 

go = uok'"^ , gi = uik~'^ . (3.4) 

Substituting the ansatz (j3.ip and evoking the Landau-gauge decouphng theorem of App. |Dl 
the FRGE (jl.ip reduces to the simple form 

OiFfc = ^g^'^" + 5^" . (3.5) 

Here, the 5S'''^^-trace captures the contributions of the physical degrees of freedom in the gravi- 
tational sector, the transverse-traceless tensor /i^j^ and the trace h, while 5™' contains all terms 
that reoccur in any action for single metric gravity, including gauge fixing, ghost terms and Ja- 
cobians for the implemented decompositions. The latter are universal, in the sense that they do 
not dependent on the gravitational part of the effective average action (j3.2|) . Their contribution 
is found in App. [E] and reads 

5""'= / d^x^ {-2Q2[f'] + Ri-l,Qi[f]-lQ2[f]+Q3[f] + 6Q,[f])} , (3.6) 
with the functions /" defined in eq. (|E.lip . 



3.1 The contribution of cS^-^^^ 

The contribution of the gravitational sector S^^^^ follows from first expanding ()3.2p to second order 
in the fluctuation field, and subsequently carrying out the TT-decomposition (j2.6p . Retaining 
the terms built from the transverse-traceless tensor /i^^^ and the scalar trace h only and dropping 
the bar from background quantities, this results in 

j.grav,quad^l J d^^^ l^h^^ [k^^^^ + M^^^'] /i^^ + [Kq + Mq] /l } , (3.7) 



with 



= -5° (niA + no) , 
Ko = fiiiA + ino , 

: 0. 



(3.8) 



Notably, there are no cross-terms and a potential is only present in the transverse part. The 
operators IK2T and M2T should be understood as restricted to the transverse-traceless subspace. 
Instead of implementing this restriction explicitly, it is more convenient to leave them in the form 
(j3.8p and employ the necessary projection operators (jB.SP inside the trace later on. It is precisely 
these projectors which lead to non-trivial terms that allow us to demonstrate the virtues of our 
algorithm. 
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3.1.1 The IR-cutoff 

We now construct the matrix-valued IR cutoff TZk for the gravitational sector, which is added to 
the action in the form 



(3.9) 



According to the Type Ic scheme (j2.9p . we only regulate the kinetic terms Kj, so that there are 
no curvature terms entering TZ^- In the case at hand, this is achieved by 



7^' 



k,2T 



n 



fe,0 



lUiRk ' 



(3.10) 



The inclusion of TZk then results in replacing 

]fC2T ^ P2T = IK2T + ^fc,2T , Kq ^ Pq = 1^0 + ^fe,0 ' (^-H) 

which, instead of A, now include the regulated Pk as arguments. 

(2) 

3.1.2 The perturbative inversion of r^^^ + TZk 

Since the truncation ansatz p.2p contains only terms to first order in the curvature, it is sufficient 
to carry out the inversion neglecting all terms of the order M^. Adapting the general formulas 
(|2.12p to independently invert the tensor and the scalar parts, we have 



n2T 



11 1 

XIoT ■ Mot ■ IIsT 

P2T P2T P2T 



and 



-1 

hh 



rifr;;;^ — Iltr • 
Phh 



(3.12) 



(3.13) 



The projection operators lis are defined in ()B.5p and ensure that all terms are restricted to their 
proper subspaces. 



3.1.3 Evaluation of the operator traces 

In the following, let us denote the trace contributions by 72, where s is the spin of the trace 
argument and i enumerates the contributions. In the scalar sector, there is only the operator 
trace, 



1 



dtRl 



while in the tensorial sector we have 



T2T = ■^Tr2 



dtRl^U2T 



P2T 



-Tr 



2T 



P2T 



dtRl^ 



'2T 



1 



1 



P2T 



U2T ■ M2T • n2T -f^ dtRl^ U2T 



1 



(3.15) 
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Here, P2T and Pq are the scalar parts of the matrix valued operators P2T and Pq 

P2T = -ui Pk - uo , Po = Pk + luo . (3.16) 

The absence of contributions from commutators involving the projection operators is ensured by 
(lAlOD . 

The scalar and the first 2T-trace are easily found with the heat-kernel coefficients calculated 
in App. IC.3I and summarized in Table [H The 7^-trace containing projector insertions is more 
complicated, but for the present truncation it is sufficient to simplify it to 



= -■;TTr2 



p2 



+ 0{R^), (3.17) 



replacing the projector with its zero-order part (jB.lSP and commuting all covariant derivatives 
freely, since M2T already contains one power of the curvature. Its explicit evaluation requires 
the off-diagonal heat-kernel scheme explained in Subsection 12.3.21 and represents a non-trivial 
demonstration of the proposed automatization technique. Using this scheme, the traces 7q^, 7^\. 
and 7^ evaluate to 

j d'xVa {Q2[f^] + lRQi[f^]} , 

j ^'^^^ {^Q^UIt] -iRQAfh]} , (3.18) 

with the dependence of the operator traces on the regularized Laplacian being captured by 

f2T ■.=^^dtRf , := j^dtR,,o ■ (3.19) 

Combining these results, the full gravitational sector contribution then reads 

+ {\Qi[fl] - lQi[fh] + fn^Q2[fh]) R 



(3.20) 



Notably, this is exactly the result obtained by earlier computations specifying the background 
metric to be the one on the four-sphere [24] . We stress, however, that in our derivation of (j3.20p 
the background metric g^y is left completely unspecified. This demonstrates the background- 
independence of the proposed algorithm and highlights its generalized applicability. 

3.2 /^-functions and non-Gaussian fixed point 

The /3-functions governing the running of Newton's constant and the cosmological constant can be 
read off by substituting ()3.20p and (13. 6p into ()3.5p and comparing the coefficients of Iq = / d'^x^ 
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and h = J d'^x^/gR. The result is conveniently written in terms of the threshold functions 

(3.21) 



T{n)J^ (z + i?(o)(z)+a;)P 

1 r...n-i R'-'X^) 



1 /■°° 

with Rkip^) = k'^ R''^\p^ / k"^) . Their relation to the Qn[f] is given by 



Qr 



Ui{Pk+Uj)P 



dt{uiRk) 



2fe2(n-P+i) (ci>P + 1 + 2)ci)^^(a;/fe2)) , (3.22 



which for the functions (j3.19p and (lE.lip yields 



Qn[/on = 2(|5i)^-^A;2("-2p+2)(^^P(2^) + i(^ + 2)$f;(|f^^ , (3.23) 

Q„[r] = 2fc2("~P+i)ci>P(0). 

Inserting these relations into (j3.20p and (j3.6p and matching the power of the background 
curvature, the RG equation leads to 

dm = - 4,0 + (^Ao(f ) + (4^ (l + i?o(f ) , 

= - 2<71 + (4V^l(f^) + (4^ (l + ) > 

which can be solved for the /3-functions for the dimensionless Newton's and cosmological constant 
{g and A) with the use of the definitions ()3.3p and (|3.4p . This yields 



dtX ={m - 2)A + 1^ (Ao(A) - iwSo(A)) 



5*5 =(r?Ar + 2)5, 
with the anomalous dimension of Newtons constant 



(3.25) 



and the cutoff shape dependent functions 

^o(A) =5$^(-2A) + <i>^(-|A) - 4$^(0) , 
Bo{\) =5$^(-2A) + $^(-|A), 

A, {\) = - l^\{-2\) + i<I>l(-|A) - ^i^l{-2\) (3.27) 

- ^$1(0) - |cD2(o) + 2$i(0) + 12c&|(0) , 

B, {X) = - |$i(-2A) + i$l(-|A) - f i'i(-2A) . 

Evaluated with the Optimized Cutoff R^^\z) = (1 — z)Q(l — z) |45j . this system has, in 
addition to the Gaussian fixed point, a non-Gaussian fixed point (NGFP), in agreement with 
previous works [35l [Ml EB [38] . This NGFP is located at 

g^ = 1.0021 , A* = 0.134414 , c/*A* = 0.134696 , (3.28) 
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and has the stabihty coefficients 

61 = 2.37141 ±i 2.27954. (3.29) 

The /3-functions (j3.25p and numerical values for the NGFP (j3.28p and (|3.29p are very close, 
but not identical to the ones obtained in [23\ 124] . This difference can be traced back to the use 
of the two different cutoff schemes (Type I vs. Type Ic), which slightly alters the contribution 
steming from the gauge-fixing part of the action. We find it very encouraging, though, that both 
cutoffs lead to the same qualitative results. 

4 Discussion and future perspectives 

In this paper, we analyzed the structure of Wetterich-Type Functional Renormalization Group 
Equations (jl.ip , which provide an exact description of the Wilsonian renormalization group flow 
of a given theory. The flow equation contains the same physical information as the path-integral 
formulation, but is generally much more manageable. The prime computational obstacle when 
extracting non-perturbative physics from this equation is the evaluation of the operator traces 
appearing on its right-hand side. 

As the central result of this paper we presented an explicit algorithm, the universal RG 
machine, that bypasses this obstacle by relating the operator traces to the off-diagonal Heat- 
Kernel expansion of a standard Laplace-Operator [29^ 130] (for reviews and further references, see 
also [311 132] )• The explicit steps in the algorithm are detailed at the beginning of Sect. [2] and 
allow for the construction of the /3-functions arsing from a given truncation by purely algebraic 
methods. Each of these steps can be handled by computer algebra software, which is highly 
convenient for the treatment of truncations beyond a certain complexity, where the /3-functions 
arise from the sum of a large number of terms. As a welcome byproduct, the construction 
lifts the technical restriction of working with a particular background, manifestly showing the 
background independence of the formalism. Thus, the algorithm provides access to information 
that is otherwise restricted by the choice of the background. 

These features open the door to a multitude of applications which up to now where out 
of computational reach. Furthermore, while our motivation for constructing the universal RG 
machine originates from studying the renormalization group flow of gravity in order to obtain 
further insights into the asymptotic safety conjecture, we stress that the applications of this 
construction are not limited to gravity and the algorithm is easily adapted to address open 
questions within other gauge theories. 

As an illustration of our algorithm, we re-derived the /3-functions of the single-metric Einstein- 
Hilbert truncation [35l [SSI [33 [31]. The resulting flow equation was evaluated for a generic 
background metric, demonstrating that the resulting /3-functions depend on the gauge-fixing and 
cutoff-scheme, but are explicitly independent of the choice of background metric. Moreover, the 
new type of cutoff, intrinsic to the universal RG machine, gives rise to the same fixed-point 
structure found within all other cutoff-schemes [38[ [2T] . 
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The construction of the /3-functions is considerably simphfied by working in the geometric 
Landau gauge ()E.3p . As demonstrated in App. [Dl this choice leads to a factorization of the 
operator traces into a gravitational sector, containing only the contribution of the transverse- 
traceless and trace parts of the metric fluctuations, and a "universal sector" which captures 
the contributions of the gauge fixing, ghost, and auxiliary terms arising from the applied field 
decompositions. Once computed to sufficient order, the latter sector can be recycled for any 
truncation of the gravitational effective average action, as it does not depend on the choice of the 
gravitational part of F^. 

A basic building block of our algorithm are the heat-kernel coefficients for transverse vectors 
and transverse-traceless matrices, which we computed on a manifold without boundary up to 
second order in the curvature. The result is given in Table [T] While independent of the choice 
of metric, the result depends on the topology of the manifold, which affects the computation 
via the zero modes of the operators implementing the transverse decomposition of a vector field. 
Based on this observation, it is natural to expect that this topology-dependence will carry over to 
the gravitational /3-functions. Thus, we expect that at a sufficient level of sophistication (which 
is beyond the Einstein-Hilbert case analyzed in Sect. [3]) the gravitational /3-functions will be 
sensitive to topological effectsif] 

Notably, the universal RG machine aids practical computations in any expansion scheme of 
the flow equation, such as a vertex expansion in a suitable background field. This has already been 
demonstrated in the context of computing the ghost-wave-function renormalization for gravity 
[49| , where the expansion is controlled by the background ghost flelds |5H HOj [50] . In a similar 
spirit, one can systematically carry out bimetric computations along the lines [13118], providing 
a handle on, e.g., the running of the gauge-parameters. 

In the case of gravity, the presented algorithm allows us to extend the explored theory space 
into hitherto uncharted terrain, facilitating the systematic derivative expansion of the gravita- 
tional effective average action |34j . This task is closely linked to determining the flxed point 
structure and number of relevant couplings of the theory. Using our results, it is now straight- 
forward to complete the study of four-derivative truncations initiated in [52l [26] , by including 
the scale-dependence of the Gauss-Bonnet term. This issue will be addressed in a forthcoming 
publication. Lastly, we note that, though requiring an improved understanding of the off-diagonal 
heat-kernel expansion going beyond the results presented here, the inclusion of the Goroff-Sagnotti 
term [53l \5^ can also be in principle implemented along the same lines. 
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A Commutator relations 

Commutator relations between covariant derivatives are an important building block of the uni- 
versal RG machine. In this section, we summarize the most important commutator relations for 
the case of gravity, which also occur in the example given in Sect. [3l 

The commutator of two covariant derivatives acting on an arbitrary tensor <^Qia2...a„ is given 

by 

n 

[D^ , Diy] 4>a-ia2...0Ln — RfMuajc'^ 'Pai...ak^ipak+i...a„ ■ (^-1) 

k=l 

Based on this, it is convenient to have commutators involving the Laplacian A = —D^D^j^ acting 
on scalars (/>, vectors (pa and symmetric matrices (pap 

[D^,A]cpp = R^^Da^p - 2R''/p Da^p - {DaR''/p)(t>p , (A.2) 

[D^,A]ct)p, = R^^DaCPpa - ^Rap\p D^<p,^p - 2 {DaR'^/^^p) (P^^p . 

For later reference we also note the double-commutators 
[[D^,A],A]0 =R^'^RjDpcP + 0{DR), 

[[D\ A] , A] <A^^ = -R^- [Ru^D^cP^p + AR^^^^,aD''cpp)p\ (A.3) 
+2R^p^'' [Ru^'DpcP^x + 4i?^(^.a£'"</'A)M] + 0{DR) . 

These basic commutator relations can be used to expand the commutator of a covariant 
derivative with an arbitrary function of the Laplacian in terms of the curvature. As a prerequisite, 
one verifies 

[Dp,e-'^] ^a,...a^=e-'^ |-s[Z)^,A] + ^[[Z)^,A],A]| cPa,...a^ + 0{DR, R^) . (A.4) 

It is straightforward to extend ()A.4p to functions of the Laplacian, 

[D^ , /(A) ] ^a,...a„ = {/'(A) [D^,A] + i/"(A) [[Dp, A] , A] + OiDR, R^) . (A.5) 

Here the prime denotes the derivative of /(x) with respect to its argument. This result can 
be established by replacing /(A) by its Laplace-transform, subsequently substituting ()A.4p and 
transforming back. 

Finally, for practical computations, it is useful to have explicit expressions for the commutators 
of the projection operators (|B.ip and (jB.Sp . For transverse vectors, one can use eqs. ()2.15p . (jA.ip . 
and (lA.Sp to show that 

[UtL ^ /(A) [UtL cP, = [UtL /(A) - /'(A) [Ct], ^ , (A.6) 



20 



where 

" = (s^m"^- + ^R'^^D^D^Dp) D" , (A.7) 
and all expressions are valid up to terms 0{R^, DR). Analogously, one can show that 

[n^T]^, /(A) [n^T],. <t>ap = [n2T]^. /(A) - /'(A) [C2t]^. <^./5 , (A.8) 

with the corresponding commutator 

[C2t],, = [Hst]^. ^"^ (-2i?(/5i'Jl)^)Z), + 4i?\/.)I),4"Z)^)) ^ . (A.9) 

In order to arrive at this expression, it is useful to observe that the bracket is generated by the 
first term in ()B.15p . All other commutators are either proportional to g^p, D^, or Djj and are 
therefore annihilated when contracted with n2T- 

A central property of the tensors ()A.7p and ()A.9P is that they are traceless with respect to 
the internal indices 

tn [CT]/ = 5i; [Ct]/ = 0, 

/ \ . (A.IO) 

tr2 [C2tI, = 1 [6i^6^, + 5^^6{) [Ct],, "^ = 0. 

This ensures that commutators of this type do not contribute to the Einstein- Hilbert truncation. 



B Transverse-traceless decomposition 

Another central ingreedient of the universal RG machine is the York-decomposition of the fluctu- 
ation fields into their transverse and longitudinal parts [551156]. In order for the decomposition to 
be well-defined, we will assume that the background is a closed (i.e., compact without boundary) 
and complete mamfoldS Under these assumptions the decomposition is umque, up to ambiguities 
associated with Killing vectors or conformal Killing vectors of the background which constitute 
zero modes of the projection operators implementing the decomposition. In the sequel, we will 
first discuss the decomposition of vector fields, before turning to symmetric tensors in Subsection 

El 



B.l Decomposition of vector fields 

Let us start with the decomposition of a generic vector field, which typically appears in the 
gravitational gauge- and ghost-sector. The split into its transversal and longitudinal part follows 
from (j2.4p and is subject to the differential constraint ()2.5p . The systematic implementation of 
this decomposition at the level of the universal RG machine requires the introduction of suitable 
projection operators on the space of (unconstrained) vectors 

[li]/ = 5/ , [Hl]/ = -D^A-^D- , [Ht]/ = [li]/ - [Hl]/ . (B.l) 

^Strictly speaking, the assumption of the background being closed can be relaxed and one could consider 
asymptotically flat backgrounds when appropriate assumptions on the fall-off of the metric and fluctuation fields 
are implemented. 
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These satisfy the natural projector properties 

nL-nL = nL, nT-nT = nT, nL-nT = nT-nL = o, (b.2) 

which also ensure the orthogonality of the decomposition. Applying these projectors to (j2.4p . it 
is then straightforward to establish that 

[nTC]^ = Cj, [U^Cl = D^r^. (B.3) 

All of these relations are exact in the sense that they do not refer to a systematic expansion 
in the background curvature. For practical computations, however, it is convenient to have an 
explicit expression for IIl, where the Laplacian is to the very left, or very right of the matrix 
structure: 

= [-A'^D^D^ + A-^R^^D^D'' + 0{R^)\ . 
These can be obtained as a perturbative series in R by making use of the commutator relations 

coil . 



B.2 Decomposition of symmetric tensor fields 

i^jy can be decomposed in a i,ia,iiovciDc-i,ia,^,cicoo pan, n-^j^. 



The gravitational fluctuations huv can be decomposed in a transverse-traceless part /ijj^, a vector 



and a scalar h encoding the trace-part. The explicit deconstruction follows from (|2.6|) with 
the component fields satisfying (j2.7p . The "standard" TT-decomposition, employed for example 
in [571 ES], and needed to eliminate all the B-terms, includes a splitting of the vector into 
its longitudinal and transverse parts, as in (j2.4p . The reason we are not implementing that here 
is that it represents a superfluous complication for the projection on the 2T-component of the 
fluctuation, and thus also for the derivation of the associated heat kernel coefficients in App. lC^ 
In order to distinguish the decomposition (j2.6p from the "standard" TT-decomposition, we will 
refer to (j2.6p as "minimal" TT-decomposition. 

When evaluating operator traces on the space of 2T-fields it is again convenient to introduce 
covariant projection operators onto the transverse-traceless, vector, and scalar subspaces. These 
can be constructed from 

V^trVr =\9<.fi9''\ (B.5) 

Here, I2 is the unit-operator on the space of symmetric matrices and 

\P^\% = 2Z)(„5^) - y^^D^ , [P^'X = i^^'u - R'. - ^D^D^Y' . (B.6) 

Applying these projectors to ()2.6p . it is easy to verify that they project h^y onto its vector- and 
trace-part respectively, 

[n2L/i]^, =[Pi]%i. = D^i, + D,i^ - Ig^^D^i^ , [HtM^, = -d9^.uh. (B.7) 
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Furthermore, they satisfy the standard projector identities 



n 



tr 1 



(Bi 



ensuring the orthogonahty of the decomposition on an arbitrary background. Based on the 
properties ()B.7p . the projector on the transverse-traceless part hj^^ is given by 



= I2 - - Htr , [Il2Th] 



(B.9) 



This last projector is a key ingredient for the construction of the heat-kernel expansion on the 
space of /ij^-fields subject to the constraints (|2.7p in Sect. IC.31 

The presence of the complicated pseudo-differential operator [^2~^]^ severely complicates 
working with the projector n2T- To bypass this obstruction, it is necessary to recast this operator 
into a power series in the background curvature. 



(B.IO) 



ip + 0{R\DR). 



In order to derive this expansion, we abbreviate q = [d — 2)/d < 1 and read the formal inverse 
()B.6j) as a geometric series in the q- and curvature terms. The perturbative inverse then assumes 
the form 

(B.ll) 



where P2,o a-nd P2,r capture the resummed contribution of the g-term at zeroth and first order 
in Ra, respectively. Explicitly, we can write 



r,— n ^ ^ 



n=0 

oo 

E 

n=0 



(B.12) 



and 

[P2,Rtip 



A"" A ' ^ A"+2 

n=l 

oo 

n=l 



1- " ^ 

A"" A 



-I^R^^^^D^DP) + "^Rf^-D^D.D^D^ 



(B.13) 



A " A^(i-1A3^(° ^4(d-l)2A4^ ^^^'^^^^ 
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Substituting them into (|B.11[) and applying the commutator (1A.5|) to the first term in (|B.12[) 
finally proves (jB.lOp . The resummations ()B.12p and (|B.13P are exact up to terms 0{R^,DR). 
Higher orders can be calculated in the same way, taking higher order commutators into account. 
With the perturbative inverse of at hand, it is now straightforward to give the projector 
in terms of a perturbative series in the background curvature. Omitting terms 0{R^,DR), 
it is convenient to express the series as 



[n2L 



1 P'^d) 



no[^ + mi, + nl 



2L 



^2L 



pu 



[nO[ + ni£ + ny^/- 



(B.14) 



^pcr J 



where, H^i^ and denote the part with no curvature terms and the Laplacians written to the 
left and right, and Ilg^ capture the corresponding commutator contributions from moving the 
Laplacians, and IlgL is the linear curvature contribution originating from the inversion formula 
(|B.10p . A straightforward computation shows 



d-2 



with the resulting commutator pieces reading 



(B.15) 



X 
d-2 



1' , 

2L J a/3 



per 



R^^Oxd"; - 2R/^^''Dx) [ I2 - ''"i, 
^ g^p R^^ D.D'^ ^ [ I2 - I, , 

^ (^Rf^a^Oxd^p) - 'iR^'(J'p)D^ + -j^^gaf}R'"'D^ 
+ (2R(^^^Dp)DxD^ - IR^^pDxD^D'' - Df^^Dp^D^R^'') 



(B.16) 



X [ I2 - n 



tr 



per 



71/ 



Finally, the commutator free piece appearing in both expressions is 



(B.17) 



Recall that, at the level of the Einstein-Hilbert truncation, all covariant derivatives in ()B.16P and 
(IB.17P can be commuted freely, as we drop all terms of order 0{R^ , DR). The expansion (IB.14P 
is the central ingredient in practical computations to handle the projection operators in a trace 
argument. 



24 



B.3 Jacobians for the decomposition of vector and symmetric tensor fields 



The decomposition of vector fields (12. 4p requires an adaption of the path integral measure by a 
Jacobian Jvec such that 

V[C^] = J^ec V[Cl]V[ri\ . (B.18) 

Following ^58j , the Jacobian can be found by explicit evaluation of a normalized Gaussian integral 
with the decomposition inserted, 



1 



^[C^] exp 



(B.19) 



P[r/] exp — ^ / d'^x^ r]Ar] 



The remaining Gaussian integral gives a factor of [det (•o)(A)]^2 for bosonic or det(o)(^) foi^ 
fermionic fields. 

Likewise, the employment of the minimal TT-decomposition (j2.6p of the metric tensor requires 
to insert an appropriate Jacobian J. It contributes a compensating factor to the functional 
integral 

-Dih^u] = Jtt V[hlMh]V[^^] . (B.20) 
This Jacobian is also found in the same way by solving Gaussian integrals. We then find 

1= J P[V] exp[-i J d'^x^h'^'h^, 
= Jtt / V[hlMh]V[i,] exp [ - i / d^x^g {h'^'^hl, + + 2i^M\e) 1 (^-21) 



cx Jtt det (i)(M) 2 , 



with 



M^', = AS'', - (1 - - R'^u . 

Solving for Jtt leads to 

Jtt ocdet(i)(M)2 =det(i)(M) 2det(i)(M)^ 

V[T^]V[b^]V[b^] exp[-i / d'^x^ (T^M^T'^ + 6^M V) 



(B.22) 



(B.23) 



where use has been made of the Faddeev-Popov trick to exponentiate M'^i, introducing new 
auxiliary bosonic (T) and fermionic (6, b) vector fields. 



C Heat-kernel coefficients for fields with differential constraints 

We now use the projection operators of App.[B]to compute the early-time heat-kernel coefficients 
for transverse vectors and transverse-traceless matrices on a generic compact manifold without 
boundaries. While the projectors are valid for any dimension d, we will restrict ourselves to the 
case d = 4: for simplicity. The generalization to arbitrary d is straightforward. 
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C.l Standard Heat-kernel coefficients for unconstrained fields 



We start this section by summarizing the standard heat-kernel expansion of the minimal second 
order differential operator A = —D^ — Q for unconstrained fields. For small values s, the early- 
time expansion becomes 

Tr^ [e"*^] = j^^^ J d'^x^/g [tr^ oq s tr^ 02 + tr^ 04 H ] , (C.l) 

with coefficients [59] 

ao =1 , a2 = P , 

Here the subscript s = 0,1,2 indicates that the operator trace Tr^ acts on the space of scalars, 
vectors and symmetric matrices, respectively, while tr^ is a trace over the internal indices. Lastly, 
P = ^Rlg + Q and TZf^i, = 2D[^D,y]. For the special case where A = {—D^ — qR)ls the traced 
heat-kernel coefficients are 



troao 


= 1, 


troa2 = 


triao 


= 4, 


tria2 = 


tr2ao 


= 10, 


tr2a2 = 



i^i?, troa4 — ^^^-^r-R^ + imiRuuaisR^'^"'^ — Rfii/R^'^) 



m+Mi? , tria4 = ^^^R^ - ^^R,.R^^ - ^R.ua^R^''"^ , (C.3) 

^^i^^i? , tr2a4 = ^^'^gg'^^ R^ — jgR^uR^^ — ^R^iuapR^'^"'^ , 

For q = 0, these are summarized in the first three columns of Table [TJ These coefficients provide 
the starting ground for computing the heat-kernel coefficients for the constrained fields in the 
following. 

C.2 Heat-kernel coefficients for transverse vector fields 

Our strategy for computing the heat-kernel for the constrained fields builds on the observation 
that it can be obtained from their unconstrained counterparts by inserting the suitable projection 
operators. For the case of transverse vectors, the corresponding projector is given by (jB.ip . 
yielding 

TriT [e-^^] = Tri [c-^^Ht] 

= TVi [e-*^] + TrJ, [D^A-^D'^e-'^] (C.4) 
= Si + So. 

Here, the prime on the scalar trace highlights that the constant scalar mode does not contribute, 
since it drops out from the decomposition (|2.4p . 

The expansion of the trace Si can be directly read off from ()C.3p . For the evaluation of 5*0, 
we first calculate the full trace, before dealing with the constant mode. The expansion of the 
trace can be found by first commuting D^^ with e~*^. Using ()A.4p together with ()A.2p and ()A.3P 
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acting on scalars and exploiting the cyclicity of the trace, this results in 

50 =Tro [Df.A-^e-'^ + sR^'^D^ + ^s'^R^^'^R^^D'')] + 0{R^) 

= - Tro [e-^^] + sTro [A-^e-^'^R'^'' D^D,] + ^Tro [A-^e-'^R'^'^R^.D^D''] (C.5) 
= To + Ti+T2 + 0{R^). 

Tq is given by the standard heat-kernel expansion of a scalar field. To evaluate Ti and T2 we 
combine the Schwinger-trick with the off-diagonal heat-kernel methods of Section 12.31 



Ti =sTro [A-^e-'^R'^'^D^D,] 

= s / dtTio e-^'+'^^R'^'D^D^ 
Jo L 



(C.6) 



lsR-^s''R^ + ^s''R^,R^^ 



~ (47rs)2 

where we substituted the i?^,y-tensor (j2.24p in the third step. Following the same route, one finds 

T2=J^ j d^x^ [-Is^R^^R^-] . (C.7) 

In order to obtain the final result, we still have to deal with the zero mode appearing in Si. Its 
subtraction from Tq is trivial, using the identity TV^ = Xro -^o"^™~'"°'^^ where Tg^e^-niode ^ _i 
a pure number, counting the constant scalar field on the manifold. In particular, it is independent 
of the volume-element. Thus, it has to contribute to an index, which in the present case is given 
by the Euler-character xe of the manifold / d^x^{R^ - m^.^R'"' + R^.ua/sR^""''^) = 
The zero-mode then gives a contribution to the heat-kernel expansion proportional to the inverse 
Euler-character, 



j,zero-mode ^ -J—-L f d'x^{R^ - AR^^R^" + R^.a^R'^'^"^) . 

(4vrj^ ZXE J 



(C.8) 



For the other terms we have to deal with undefined expressions (zero over zero), so we need to 
regularize them. We take such a regularization to be a sequence (continuous or discrete) of scalar 
modes (/>(™') such that lirrim^o't'^"^^ = 4'^^'^ = const., and such that for any m > there exists a 



well defined tensor s|j^^ for which 



A-iZ?^Z?,0(™) = 4^)</.("^) . (C.9) 



By counting the number of derivatives we assume that for a suitable regularization the limit 
sfl} = limm^o will exist, and it should satisfy DpS'^y = 0, together with g^^ s"p} = — 1. There 
is only one tensor satisfying these requirements, namely 

s!S = -^ffM.- (c.io) 

Comparing these contributions to (jCSh we observe that they are of higher order in the curvature, 
and thus cannot contribute to the heat-coefficients computed here. 
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Combining all these partial results via (1C.4P yields the heat-kernel coefficients for the trans- 
verse vector-fields 



(C.ll) 



triTao=3, triTa2 = 3^, 
This proves the fourth column of Table [TJ 

C.3 Heat-kernel coefficients for transverse-traceless matrices 

The computation of the heat-kernel coefficients in the 2T-case is completely analogous to the 
IT-case, but technically more involved. Using the projectors (|B.9p and (jB.Sp we find 



(C.12) 



Tr2T [e-^^] = TY2 [e-^^nsT] =1^2 [e"^^] - TV2 [e'^^HsL] - TV2 [e'^^Htr] 

= S2 + Si + Sq . 

Noting that 

5o = -TV2 [e-^^Htr] = -TVo [e"^^] , (C.13) 

both 5*2 and 5*0 reduce to standard heat-kernel formulas and their contribution can be read off 
from ([031) . 

In order to obtain 5*1, we first substitute n2L, (jB.SP and use the cyclicity of the trace to write 



Si = Tr 



Z)^)[l2-ntr],^/^^e-^ [Pi];, [P^^t 

11/3' 



7/3 ^ IJ/x!/ 

Tri [e-^^] +Tri [-se-'''[D\D''] [Pi]^^ [P^ 

^e-^^ [[D\D\D^] [Pi]:;,[P2-^]^ 

To + Ti+T2 + 0{R^). 



+ Tri 



(C.14) 

+ 0{R^) 



The commutator in the second step follows from (IA.4P and we exploited the orthogonality of the 
projectors, which ensures Iltr-Pi = 0. The commutators ()A.2p . (jA.Sp and the perturbative inverse 
P^^, (lBlO|) . can then be explicitly substituted. Expanding the products, the full trace is then 
built from linear combinations of basis integrals, which can be evaluated with the off-diagonal 
heat-kernel along the lines of (|C.6p . For convenience, these are summarized in Table [2j 

At this point, the following technical remark is in order. When evaluating operator traces 
with the off-diagonal heat-kernel it is crucial that all covariant derivatives appear in a totally 
symmetric way. This requires their symmetrization via D^Di, = D(^^D^^ + ^[D^, D^], with the 
antisymmetric part expressed as curvature tensors. In this respect it is important to note that the 
open vector indices /3 and /3' appearing in ()C.14p may not be contracted immediately. Instead, f3' 
must be thought of as contracted with a vector c/)^/ to the right of the expression. Only once all 
commutators have been expressed through curvature tensors, the "open indices" f3 and /3' may 
be contracted. 
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sTrn 




sA^21 






LA 


sA^/.!/^ n 1 


= 1^ f d^x^/q \ } R }^R^ + iRuvR^""] 

(47r)2 J -^Vy L 4s 12 ' 6 A"^ J 


sTro 








sTVo 












-^^R'^'^R'^^D^D^D^Dp] 




sTro 




-^^R^pR^'^'^^D^D,] 





Table 2: Typical operator traces arising in the computation of perturbed traces. These are 
evaluated using the off-diagonal heat-kernel methods discussed in Subsection 12.31 



Following this route, a lengthy but straightforward computation yields 



^sR + ^s'^R^ - '^s^R^.R^^- - R pM-P- 
3 27 27 ^ 3 



and 

^17 



-.(-,'5 R^yR^ ~^ o RfiupuR^ ^ 
io O 



(C.15) 



(C.16) 



(Atts) 

Combining this result with the standard vector-trace Tq, we arrive at the vector contribution to 
the 2T-trace, 

Si = I (fx^ [4 + IsR + §^2^2 ^ ^^2^^^^^. _ 131^2^^^^^^^.p.] _ 

As observed above, 5*2 and are standard heat-traces involving symmetric matrices and 
scalars, respectively, so that they can directly be inferred from ()C.3p . Adding ()C.17p . the heat- 
kernel coefficients for the 2T-trace are obtained as 

tr2Tao = 5, tr2Ta2 = -|i?, iT2TaA = -^R^ - ^R^^.R^" + ^R^^paR^'P'' . (C.18) 

This constitutes the main result of this subsection and proves the fifth column of Table [TJ 

We close this subsection with a remark on the zero modes, which do not contribute to the 
minimal TT-decomposition and have to be subtracted from the coefficients. Inspecting (j2.6p . we 
conclude that vector fields satisfying the conformal Killing equation 

D^iu + D,i^ - = (C.19) 

will not contribute to the fluctuation field, and therefore have to be excluded from the operator 
traces. The number of modes to be excluded from the trace is then given by the number of Killing 
vectors tt-kv (satisfying D^£^y + Dy£^^ = 0) plus the number of conformal Killing vectors ncKV 
which solve (|C.19P with D°^^^^ ^ 0. Following the logic for the zero modes of the previous 
subsection, these modes give rise to a correction of Si proportional to the Euler integrand 

traxaf "^'"^^"^^ = ^ (r^ - AR.^R'^'' + R^.^^R^'''^^) , (C.20) 
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where = raxv+'^CKV is the number of independent solutions of (jC.19|) . This extra contribution 
is absent in the generic case where the manifold does not posses particular symmetries. However, 
it plays an important role when comparing our general computation to the spherical result in the 
next subsection. 

C.4 Comparison of the Heat-kernel coefficients on a spherical background 

An important check for the heat-kernel coefficients (jC.lip and (|C.18P is provided by the special 
case where gfj_y is the metric of the four-sphere S^. In this case, the coefficients shown in Table [T] 
have to reproduce the earlier results [MIEH]- In this subsection, we verify that this is indeed the 
case. 

Since is a maximally symmetric space, all its curvature tensors can be expressed in terms 
of the metric and the covariantly constant Ricci-scalar 

1 1 

Rp,v = ^ RdiMUi Rfivpa ~ Y2 ^ {dppQucF ~ QpaQvp)- (C.21) 

These relations imply that the curvature-square and volume-terms can be expressed in terms of 
the Ricci-scalar, 



R^.R^"' ='^R^ , R^.p^R'^'P'' = ^R\ V0I54 = 3847r2i?-2. (c.22) 



Moreover, the four-sphere has xe = 2 and admits nxv = 10 Killing- and ncKV = 5 conformal 
Killing vectors, totallying to = 15. Substituting this data into (jClip and ()C.18p . taking into 
account the contributions arising on a special background with symmetries ()C.20p . one obtains 

triTaol54=3, triTa2|s4 = j-R, tiiTailgi = -j^R"^ , 

5 1 2 (C.23) 

tr2Tfflol54 = 5, tr2Tfl2|54 = —qR, tr2Tfl4|54 = "432-^ • 

These are exactly the tabulated results given in [361 138]. We thus establish that the projector- 
method employed here recovers the heat-kernel results obtained previously. 

D Decoupling gauge degrees of freedom, ghosts, and auxiliaries 

The use of standard gauge fixing terms and field decompositions leads to recurring terms in 
the RG equation that can be evaluated independently of the specific model under consideration. 
Here, we will demonstrate that it is possible to separate the contributions of the gauge fixing 
and Faddeev-Popov terms to the running of coupling constants from those of the physical part 
of the action. To achieve this, a gauge field A is first decomposed in background and quantum 
fluctuations as A = A + a, and the quantum fields are then split into component fields such that 
the gauge- fixing condition contains only certain modes, denoted as Cg. With an appropriately 
chosen gauge fixing condition F(A, a) = F{A, Og) linear in the fluctuation field, the gauge fixing 
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term can be written in the form 



(D.l) 



for some operator G depending on the background fields only. It should be noted that invariance 
with respect to the background gauge transformation can always be maintained irrespective of 
the separation of modes. 

In order to write the structure of the RG equation in terms of the decomposed gauge fields, 
we can arrange the latter in a multiplet <I> = (i?i>, Og), where (j) represents all remaining component 
fields of A and other fields occurring in the action like the Faddeev-Popov ghosts. Any of these 
fields may be later further decomposed into component fields. 

With the gauge part formally separated, the quadratic part of the average effective action 
assumes the form 



nquad 



(D.2) 



2 [q IG + H 

where the block structure is chosen such that the lower right is the quadratic part of the pure 
gauge component, L denotes the contributions from all remaining fields, and Q and Q are mixed 
terms. Schematically, the FRGE (jl.ip then becomes 

-1 



dtTk 



1 



STr 



dt 



Rl Rq 
R^ IRg + Rh 



(D.3) 



L + Rl Q + Rq 

Q + Rq ^{G + Rg) + H + Rh 

The crucial thing to note here is that a suitably constructed cutoff for the gauge component Rq 
according to the scheme (j2.8p must come with the same a-dependence as G. In (jP.Sp . this feature 
is highlighted by displaying the a-dependence explicitly. A physically motivated choice of the 
gauge parameter is that of the Landau-de Witt gauge (a = 0) , since there is no smearing of the 
gauge condition. Moreover, a = is expected to be a fixed point for the gauge parameter [601I6T] . 
Thus a = is a natural choice. The occurrence of in the cutoff term does not allow to take 
the limit a — )• right away, however. Rather, we must keep track of a in linear order to cancel 
the inverse. 

With the help of the general expression for the inverse of a block matrix (j2.1ip . we can expand 
the inversion around a = up to linear order, yielding 



L Q 
Q tG + H 



- -1 


o" 








+ a 










'QG-^ 



(D.4) 



L-^QG-^QL-^ 
-G-^QL-^ 

This formula requires the invertability of L and G, a fact which we have to keep in mind in the 
explicit construction of the gauge fixing term. 

Multiplying ()D.4p with the cutoff matrix and taking the limit a — >• 0, we find 

1^^ r {L + RL)-^dtRL 
2 [ {G + RG)-^dtRG 

4sTr [(L + RL)-^dtRL\ + ^STt [{G + RG)-'dtRG] ■ 



ar 



k 



+ Oia) 



(D.5) 
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In the Landau gauge limit, the structure of the RG equation simplifies such that the mixed terms 
Q and Q as well as the block H drop out identically. Therefore, the pure gauge components 
decouple from the remaining fields and their contribution is determined by the gauge fixing term 
only. This feature will be exploited when computing the universal contribution to the gravitational 
RG flow in the next section. 

Furthermore, a standard ghost term of the form 

rf = j d'^x./^C^-^5cA = j (fxy/^CMC (D.6) 

is needed in the case of non-Abelian gauge theories. Because of its simple quadratic form, it is 
possible to examine the structure of L by separating the ghost fields from the remaining ones. 
Making use of the background independence of the RG equation, we can set the background ghost 
field to zero, eliminating the cross-terms stemming from the ghost-gauge field interaction. Note 
that this choice will of course not allow us to extract the running of the coupling constants in the 
ghost sector. The block structure then becomes 

K {) " 

L= -M , (D.7) 
M _ 

where M is the kernel of the ghost action ()D.6P and K is the still undetermined block containing 
all other second variations of the action. Since the cutoff Rl must be of the same form, the 
resulting trace in the RG equation decomposes as 

^STr(L + RL)-'dtRL = ^Tv{K + RKy^dtRx - Tr(M + RM)-^dtRM , (D.8) 

and thus 

dtTk = ^Tr{K + RK)-^dtRK + ^Tr(G + Rar^dtRG - Tr(M + RM)-^dtRM ■ (D.9) 

In order to arrive at an entirely self-contained treatment of those degrees of freedom, it 
is necessary to include a sector of auxiliary fields that formally correct the counting of zero 
modes of the applied decompositions. Using the Faddeev-Popov trick, these terms are found 
by exponentiating the Jacobian created by the corresponding projection operators. Since these 
operators always act linearly, the resulting terms will be strictly quadratical in the auxiliary fields 
and decouple from the other fields like the ghost sector by choosing a zero background. 

Therefore, the gauge and ghost as well as the auxiliary terms decouple from the physical field 
components in the Landau-de Witt gauge. This result allows us to evaluate their contribution to 
the RG flow completely independently, leaving the matter fleld content of the model arbitrary, 
and to give universal expressions to be reused in a broad class of models. 

E Universal contributions to the gravitational RG flow 

In this appendix, we compute the re-occurring contributions S^^^, arising from the gauge fixing, 
ghost, and auxiliary actions, up to linear order in the background curvature. These naturally 
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organize themselves according to 

^uni ^ ^gf ^ ^gh ^ ^aux ^ 

and we will evaluate these traces in turn. Again, we set d = 4 for simplicity, but the generalization 
to arbitrary d is straightforward. 

E.l Gauge fixing the gravitational action 

The construction of the effective average action employs the background field method, which 
breaks the quantum gauge transformations, but retains background gauge-transformations as an 
auxiliary symmetry (see [12] for more details). Neglecting the fc-dependence, this can be achieved 
by adding the classical gauge-fixing term 

= ^ y d'^xVd r^F^F, , = D'^h^, - l3D^h (E.2) 

to the gravitational action. Notably, S^^ is quadratic in the fluctuation field h^y. Substituting 
the minimal-TT-decomposition (j2.6p . we find 



F^ = D% + iD^D'^C. + + (i - P)Df,h , (E.3) 

where /? is a freely adjustable gauge parameter. To ensure the invertibility of the Hessian from 
the gauge fixing sector, /3 = 1/4 has to be chosen. In this case, all gauge degrees of freedom are 
captured by the vector By virtue of the decoupling theorem of App. [Dl this field decouples 
from the gravitational part of the action once the Landau limit a — )• is invoked. 

Setting /? = 1/4 and dropping the bar from background quantities to simplify the notation, 
eq. ()E.2p becomes 

S^' ~ J d^x^i^ [(^'^^A - ID^'D^ - i?^,] [,^-,A - \D^D, - T • (E.4) 

Following Sect. [21 we eliminate the differential operators of non-minimal form by the transverse 
decomposition of the vector C/^ = '^J + F)^(7 with -C^Cj = 0. In terms of these component fields 
the action becomes 



(E.5) 



with kinetic terms 



[fy^ = A25^, [Kf]=|A3, (E.6) 



and vertices 



[M^^Y = -'iR^''Dy/\-2/\Ri"'Dy + 2Ri"'Ry°'Da 
[Mf\l = 3AD^'R^, + 2D>'R^,A-2D>'R^o.R° 



(E.7) 
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(2) 

Hence, the gauge fixing sector of assumes a two-by-two block structure whose entries do not 
contain B-terms. 

Following the rule (j2.9p of regularizing the kinetic terms, the cutoff-operators become 



7^■ 



;f,lT 



gf .0 _ 9(r)3 



Pi - A' 



(Ei 



The corresponding regularized kinetic terms are denoted by Pf^ = Kf^ + TZ^'^ , where i = 0, IT. 

The final step consists of evaluating the resulting operator traces along the lines of Sect. 
[2j Since Mx and are both linear in R, the cross-terms do not contribute to the present 
truncation. Terminating the perturbative inversion at linear order in M, the IT-sector evaluates 
to 



■ IT 



Tri 



(47r)2 

while the scalar sector gives 



gf IT J- gf £ k 
IT ^IT 

d^x^ [3Q2[/']+i?QQi[/'] + 3Q3[/'])] , 



5!'=^Tro 


pgi " 


-^Tro 




. 





igf ^ pgf 



1 



(47r)2 

Here the argument of the Q-functionals is defined as 



r 



{Pk 



zdtRk 



The complete contribution of the gauge-sector is then given by 



1 



d^x^ [| Q2[/'] + R {kQi[f^] + 3Q3[/'] + 6Q4[/'])] 



(47r)2 

E.2 The ghost sector 

The ghost term corresponding to the gauge fixing ()E.3p with /3 = 1/4 is given by 



5. 



gh 



d'x^gC.r''^^ 



with the Faddeev-Popov determinant 

Its part quadratic in the fluctuation flelds is 



s: 



quad 



(E.9) 



(E.IO) 



(E.U) 



(E.12) 



(E.13) 



(E.U) 



(E.15) 
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where we again identified g = g. Since the operator (jE.lSP still contains B-type terms, we carry 
out the transverse decomposition of the anti-ghost field = + D^f) and likewise for the 
ghost. The resulting S^^^ assumes a two-by-two block form without B-type operators, 



^quad 



(fx 



C. 



with kinetic terms 



and vertices 



:f !^ + Mf !^ 



itJ ^ — ^ ' 



+ 



rj + rj 



3^2 
2'^ 



If 



if 



(E.16) 
(E.17) 



4]: 



if 



If 



2D^'R 



(E.18) 



The IR-cutoff following from is 



7^' 



gh,0 



'Pi - A^) 



(E.19) 



The evaluation of the resulting operator traces then proceeds completely analogous to the 
previous subsection. The complete ghost-trace decomposes into its IT-part 



-gh 

'it 



Tr 



IT 



1 



IT 



Tri 



1 



(FIt)^ 



_ _ 1 

and the scalar contribution 



d''x^{3Q2[f'] + R ilQiif] + lQ2[f]) } , 



(E.20) 



Trn 



1 



3>gh 




gh,0 
k 



Trr 



-4— 5^7^f '° Mf 



d'x^ { 2Q2[f] + R (iQi[/i] + lQ3[f]) } 



(E.21) 



Here, the minus stems from the supertrace for ghost-fields, and the functions are given in 
(jE.lip . The full contribution of the ghost-sector results from combining ()E.20p and ()E.2ip 



^gh 



(4^ 



d'x^{5Q2[f] + R {^Qi[f] + §Q3[f]) } 



(E.22) 



E.3 Jacobians of the TT-decomposition 

The final contribution to the universal sector comes from the Jacobians arising from the transverse 
decomposition of the fluctuation flelds. These occur in two types, the Jacobian from the minimal- 
TT-decomposition and from a vector T-decomposition, S^^^ = S^r^ -\- S^^^. The flrst one arises 
from the path-integral over the auxiliary fields given in ()B.23p . Notably, the operator coincides 
with the one arising in the ghost sector. Thus 



caux 1 cgh 

'-'mTT — 2'-' ' 



(E.23) 
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where S^^ is given by (ID.Gp . and the prefactor takes the combination of bosonic and fermionic 
vector field into account. 

Secondly, we have split two bosonic vectors and two complex fermionic vectors into their trans- 
verse and longitudinal components, which in total result in the Jacobian Jvec = [det (•o)(A)]^^. 
Following (]B.19p . this determinant can be exponentiated by introducing a complex scalar field 
s, s 

J^g^ = j V[s] exp [-sAs] . (E.24) 

The corresponding contribution to the RG flow is given by a standard scalar trace and easily 
evaluated to 

= I d'x^ {Q2[f] + iRQiif]} . (E.25) 
The trace 5"" is now found by adding (lKT2l) . (1^221) . ^K23^ . and (lOSll and is given in dMl). 
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